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Abstract. We report on the microwave surface resistance of two polycrystalline Mg 11 B2 samples; one 
consists of pristine material, the other has been irradiated at very high neutron fluence. It has already 
been reported that in the strongly irradiated sample the two gaps merge into a single value. The mw 
surface resistance has been measured in the linear regime as a function of the temperature and the DC 
magnetic field, at increasing and decreasing fields. The results obtained in the strongly irradiated sample 
are quite well justified in the framework of a generalized Coffey and Clem model, in which we take into 
account the field distribution inside the sample due to the critical state. The results obtained in the pristine 
sample show several anomalies, especially at low temperatures, which cannot be justified in the framework 
of standard models for the fluxon dynamics. Only at temperatures near T c and for magnetic fields greater 
than 0.5H C 2(T) the experimental data can quantitatively be accounted for by the Coffey and Clem model, 
provided that the upper-critical-field anisotropy is taken into due account. 

PACS. 74. 25. Ha Magnetic properties - 74.25.Nf Response to electromagnetic fields (nuclear magnetic 
resonance, surface impedance, etc.) - 74. 25. Op Mixed states, critical fields, and surface sheaths 



1 Introduction 

It has widely been shown that the superconducting prop- 
erties of MgB2 are strongly related to the two-gap struc- 
ture of its electronic states [TP ,3,4,5,6]. The smaller su- 
perconducting gap, A„, arises from the quite isotropic 7r 
bands and the larger one, A a , from the strongly anisotropic 
a bands [7]. Due to the different parity of the a and 
7T bands, inter-band scattering of quasiparticles is much 
smaller than the intra-band one, making the two super- 
conducting gaps quite different though they close at the 
same T c 0. 

According to the theory of multi-band superconduc- 
tivity [7] , the inclusion of defects would increase the inter- 
band scattering and, consequently, change the relative mag- 
nitude of the different gaps. In order to carry out inves- 
tigation on this topic, essentially two methods have been 
used to insert defects and/or disorder in MgB2: chemical 
substitution and damage by irradiation p [T0lfTTlfr2l[T3] . 
In any cases, the inclusion of defects, besides to change 
the inter-band scattering, might increase the upper crit- 
ical field and the critical current, strongly affecting the 
fluxon dynamics. Very recently, the effects of neutron ir- 
radiation have extensively been investigated on polycrys- 
talline Mg n B 2 plfT4p5|[T6l[T7]. It has been shown that 
irradiation leads to an improvement in both critical field 
and critical current density for an exposure level in the 
range 1 -f- 2 x 10 18 cm~ 2 . On further increasing the neu- 



tron fluence, all the superconducting properties, such as 
T c , H C 2, J c , are strongly suppressed. Furthermore, mea- 
surements of specific heat, as well as point-contact spec- 
troscopy, have shown that in the sample irradiated at the 
highest fluence (1.4 x 10 20 cm~ 2 ) the two gaps merge into 



a single value [i5|[T6]. 

Despite the large amount of experimental and theo- 
retical work done on MgB2, some arguments are still un- 
der discussion, such as the effects of the magnetic field 
on the superconducting properties [!8llT9ll20] . The main 
difficulty to quantitatively discuss the mixed-state prop- 
erties of MgB2 arises from the unusual flux-line properties 
due to the different coherence lengths, £ CT and S;^, asso- 
ciated with A a and A v |21j . Scanning tunnelling spec- 
troscopy on MgB 2 single crystals along the c-axis, which 
probes mainly the it band, has highlighted a core size much 
larger than the estimates based on the measured _ff c2 val- 
ues, as well as a significant core overlap at fields much 
lower than the macroscopic H C 2 [22j . Furthermore, mea- 
surements of neutron scattering from the vortex lattice 
have highlighted a spatial rotation of the vortex lattice for 
applied magnetic fields in the range 0.5-^1 T 23 . Accord- 
ing to point-contact-spectroscopy experiments |24 1I251[2"6"] , 
these unusual properties have been ascribed to the strong 
suppression of the superconductivity in the ir band, oc- 
curring in that field range. At low fields, each vortex has 
a composite structure, with er-band quasiparticles local- 
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ized in a region of radius £ CT and 7r-band quasiparticles in 
a wider region of radius In the field range 0.54- 1 T (at 
low J 1 ), the giant cores start to overlap; when the magnetic 
field is large enough to suppress the 7r-band gap, the vortex 
cores shrink and the 7r-band quasiparticles are widespread 
in the whole sample [T9"ll22| . This field- induced evolution 
of the vortex lattice is expected to affect the vortex- vortex 
and vortex-pinning interactions, making the description of 
the properties involving the presence and motion of flux- 
ons very difficult. 

The investigation of the microwave (mw) surface imped- 
ance, Z s = R s + iX s , in superconductors is a useful tool 
for determining several properties of the superconducting 
state. In the absence of static magnetic fields, the variation 
with the temperature of the condensed-fluid density de- 
termines the temperature dependence of Z s . On the other 
hand, the field dependence of R s in superconductors in 
the mixed state is determined by the presence of fluxons, 
which bring along normal fluid in their cores, as well as the 
fluxon motion [27ir25ra[50l [5T | [52] . So, investigation of the 
magnetic-ficld-induced variations of the surface resistance 
provides important information on the fluxon dynamics. 

Several studies of the mw response of MgB2 reported 
in the literature have shown that the experimental results 
cannot be accounted for in the framework of standard 
theories [33,34,35,36,37,38,39]. The temperature depen- 
dence of the mw conductivity, at zero DC magnetic field, 
has been justified considering the coexistence of two dif- 
ferent superconducting fluids, one related to carriers living 
on the a band and the other to carriers living on the 7r 
band [33|f34] . The magnetic field dependence of the surface 
resistance has shown an anomalous behavior especially at 
low temperatures; several authors have highlighted un- 
usually enhanced field-induced mw losses at applied mag- 
netic fields much lower than the upper critical field 35, 
36,37,38^. Sarti et al. [39j, investigating the mw surface 
impedance of MgB2 film, have shown that at low fields, 
when the contribution of the 7r-band superfluid cannot be 
neglected, the magnetic-field dependence of the real and 
imaginary components of the surface impedance exhibits 
several anomalies. Furthermore, a magnetic hysteresis of 
unconventional shape has been detected in the R S (H) 
curves [3"8"ll4Tj] . All these results have suggested that in a 
wide magnetic-field range the standard models for fluxon 
dynamics fail when applied to MgB2. 

In this paper, we report on the microwave surface resis- 
tance of two of the polycrystalline Mg 11 B2 samples stud- 
ied in Refs. [B l lH | ll3 j il^ | irf| . We have investigated the 
unirradiated sample, which clearly shows two-gap super- 
conductivity, and the sample irradiated at the highest neu- 
tron fluence, in which the two gaps merge into a single 
value. The investigation has been carried out with the 
aim to compare the results obtained in two-gap and one- 
gap MgB2 superconductors. To our knowledge, the mw 
response of neutron-irradiated MgB2 samples has not yet 
been investigated. The mw surface resistance has been 
measured as a function of the temperature, in the range 
2.5 4 40 K, and the DC magnetic field, from to 1 T, at 
increasing and decreasing values. We show that the results 



obtained in the strongly irradiated sample can quite well 
be justified in the framework of standard models, in the 
whole ranges of temperatures and magnetic fields investi- 
gated. On the contrary, the results obtained in the pris- 
tine sample cannot thoroughly be justified. In particular, 
the R S (T) behavior at zero field has been accounted for, 
in the framework of the two-fluid model, assuming a lin- 
ear temperature dependence of the normal and condensed 
fluid densities. At low temperatures, the field dependence 
of R s has shown several anomalies, among which a mag- 
netic hysteresis having a unexpected shape. At temper- 
atures near T c and applied magnetic fields greater than 
f» 0.5iJ C 2(T), the results are well accounted for in the 
framework of the Coffey and Clem model [30] , with flux- 
ons moving in the flux-flow regime, taking into account 
the anisotropy of the upper critical field. 



2 Experimental apparatus and samples 

The microwave surface resistance, i? s , has been investi- 
gated in two bulk MgB 2 samples. The procedure for the 
preparation and irradiation of the samples is reported 
in detail elsewhere [13,16 . The samples have been pre- 
pared by direct synthesis from Mg (99.999% purity) and 
crystalline isotopically enriched n B (99.95% purity), with 
a residual 10 B concentration lower than 0.5%. The use 
of isotopically enriched U B makes the penetration depth 
of the thermal neutrons greater than the sample thick- 
ness, guarantying the irradiation effect almost homoge- 
neous over the sample. Several superconducting properties 
of the samples have been reported in Refs. (T3UH11T5, 16, 
If 7j . For simplicity and ease of comparison, we label the 
two samples as in Ref. [T3], i.e. P0 (pristine Mg 11 B 2 ) and 
P6 (irradiated at the highest neutron fluence). According 
to point-contact spectroscopy |I5j and specific-heat mea- 
surements [TBJ, sample P0 shows a clear two-gap super- 
conductivity; in sample P6 the irradiation process at very 
high fluence (1.4 x 10 20 cm -2 ) determined a merging of 
the two gaps into a single value. 

Sample P0 has a nearly parallelepiped shape with w rs 
3.1 mm, t « 1.5 mm and h rs 3.2 mm; it undergoes a nar- 
row superconducting transition with T° nset rs 39.0 K and 
AT C rs 0.2 K (from 90% to 10% of the normal-state resis- 
tivity); its residual normal-state resistivity is p(40 K) ri 
1.6 fifl cm and the residual resistivity ratio RRR « 11, 
the critical current density at zero magnetic field is J c q ri 
4 x 10 5 A/cm 2 , and fi H c2 (5 K) rs 15 T; the anisotropy 
factor of the upper critical field at T = 5 K is 7 rs 4.4 
EH]. 

Sample P6 has a nearly parallelepiped shape with w ~ 
1.1 mm, t rs 0.8 mm and h rs 1.4 mm. The main char- 
acteristic parameters of sample P6 are: T° nset rs 9.1 K, 
AT C rs 0.3 K, RRR rs 1.1, p(40 K) rs 130 ^ cm. The crit- 
ical current density at T — 5 K and at zero magnetic field 
is J c o ~ 3 x 10 4 A/cm 2 ; it exhibits a monotonic decrease 
with the magnetic field, following roughly an exponential 
law. The upper critical field is isotropic and its value at 
T = 5 K is hqH c2 ~ 2 T. 
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The effects of the neutron irradiation on both super- 
conducting and normal-state properties of a large series 
of Mg 11 B2 bulk samples, including sample PO and P6, 
have extensively been investigated in Refs. [T3 l fT4 l [T7] . On 
increasing the neutron fluence, it has been observed a 
monotonic decrease of T c and an increase of the residual 
normal-state resistivity p{T c ). Nevertheless, it has been 
shown that the irradiation does not affect the variation 
of the normal-state resistivity Ap — p(300 K) — p{T c ). 
As suggested by Rowell [41], just Ap is a parameter that 
gives information on the grain connectivity. The results 
reported in Ref. [13] show that Ap remains of the order 
of 10 pfl cm over the whole range of irradiation level; in 
particular, in sample PO Ap = 16 pQ cm and in sample 
P6 Ap = 12 pQ cm, indicating that thermal-neutron irra- 
diation does not affect the grain-boundary properties. On 
the other hand, it has been shown that either neutron ir- 
radiation or He-ion irradiation [11] do not affect the grain 
connectivity, even at high irradiation levels, contrary to 
what occurs using heavy-ion irradiation [12J . Recent stud- 
ies by transmission electron microscopy have highlighted 
that neutron irradiation in these samples creates nano- 
metric amorphous regions (mean diameter ~ 4 nm) in 
the crystal lattice, whose density scales with the neutron 
dose ;17 . Studies on the field dependence of the critical 
current density have shown that at moderate neutron- 
flucnce levels (< 10 19 cm -2 ) such defects introduce new 
pinning centers, leading to an improvement of the criti- 
cal current density; on the contrary, for neutron fluence 
higher than 10 19 cm -2 (as for sample P6) these nanomet- 
ric defects do not act as pinning centers because they are 
smaller than the coherence length [IT]. Moreover, these 
studies have shown that in the pristine and the heavily ir- 
radiated samples the pinning mechanism is ruled by grain 
boundaries. The defects induced by neutron irradiation 
act as inter- and intra-band scattering centers; the intra- 
band scattering causes a reduction of the electron mean 
free path and is responsible for the growth of the normal- 
state resistivity. The reduction of T c has been ascribed 
to both the scattering processes and the smearing of the 
electron density of states near the Fermi surface [Tl"lll3| . 

Although in the two samples AT C is roughly the same, 
in sample P6, due to the reduced T c value, AT C /T C ms 
0.03, affecting noticeably the temperature dependence of 
the mw surface resistance near T c . On the other hand, 
from AC susceptibility measurements at 100 kHz, we have 
found that the first derivative of the real part of the AC 
susceptibility can be described by a Gaussian distribution 
function of T c , centered at T c q — 8.5 ±0.1 K with ctt c = 
0.2 ±0.05 K. In the following, we will use this distribution 
function to quantitatively discuss the results obtained in 
sample P6. On the contrary, for sample P0, AT C /T C is one 
order of magnitude smaller, not noticeably affecting the 
R S {T) curve. 

The mw surface resistance has been measured using 
the cavity-perturbation technique 27J. A copper cavity, 
of cylindrical shape with golden-plated walls, is tuned in 
the TEon mode resonating at oj/2it as 9.6 GHz. The sam- 
ple is located in the center of the cavity, by a sapphire 



rod, where the mw magnetic field is maximum. The cav- 
ity is placed between the poles of an electromagnet which 
generates DC magnetic fields up to poH «1T. Two ad- 
ditional coils, independently fed, allow compensating the 
residual field and working at low magnetic fields. A liquid- 
helium cryostat and a temperature controller allow work- 
ing either at fixed temperatures or at temperature varying 
with a constant rate. The sample and the field geometries 
are shown in Fig. la; the DC magnetic field is perpendic- 
ular to the mw magnetic field, J? w . When the sample is in 
the mixed state, the induced mw current causes a tilt mo- 
tion of the whole vortex lattice [31] ; Fig. [TJd schematically 
shows the motion of a flux line. 



(a) 



. Ha 



Ho 
-> 



(b) 



t 

Fig. 1. (a) Field and current geometry at the sample surface, 
(b) Schematic representation of the motion of a flux line. 



The surface resistance of the sample is given by 



R s = r 



i 



1 

Qu 



where Ql is the quality factor of the cavity loaded with the 
sample, Qjj that of the empty cavity and r the geometry 
factor of the sample. 

The quality factor of the cavity has been measured by an 
hp-8719D Network Analyzer. The surface resistance has 
been measured as a function of the temperature, at fixed 
values of the DC magnetic field, and as a function of the 
field, at fixed temperatures. All the measurements have 
been performed at very low input power; the estimated 
amplitude of the mw magnetic field in the region in which 
the sample is located is of the order of 0.1 pT. 



3 Experimental results 

Figure [2] shows the temperature dependence of the sur- 
face resistance in the pristine (a) and irradiated (b) MgB2 
samples, at different values of the DC magnetic field. In 
order to disregard the geometry factor, and compare the 
results in samples of different dimensions, we have normal- 
ized the data to the value of the surface resistance in the 
normal state, R n , at T — T° nset . The results have been 
obtained according to the following procedure: the sam- 
ple was zero-field cooled (ZFC) down to low temperature, 
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then Hq was set at a given value and kept constant during 
the time the measurement has been performed. 

On increasing Hq, the R S (T) curves broaden and shift 
towards lower temperatures; however, the effects of the 
applied magnetic field is different in the two samples. Al- 
though the value of H c i of sample PO at low tempera- 
tures is one order of magnitude larger than that of P6, 
the field-induced variations of R s in the two samples have 
roughly the same magnitude. In sample PO, one can ob- 
serve an anomalously enhanced field-induced broadening 
of the R S (T) curve, which extends down to the lowest tem- 
perature. On the contrary, in sample P6 the larger shift 
of the R s (T) curve induced by H is expected because of 
the lower H C 2 value. 

■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

r o m H=0 T Sample PO - 




Q-uJ I I 1 

5 10 15 20 25 30 35 40 45 

I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I I I I | I 

~ o n H=Q T Sample P6 t 
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23456789 10 



T(K) 

Fig. 2. Normalized values of the surface resistance as a func- 
tion of the temperature, obtained in the two samples, at dif- 
ferent values of the DC field. R n is the surface resistance at 
T = T C . 

The field-induced variations of R s have been investi- 
gated for different values of the temperature. For each 
measurement, the sample was ZFC down to the desired 
temperature; the DC magnetic field was increased up to 
a certain value and, successively, decreased down to zero. 
Figures [3] [4] and [5] show the field-induced variations of R s 
for the two samples, at different temperatures. In all the 
figures, AR s {Hq) = R S (H ,T) — R res , where R res is the 
residual mw surface resistance at T = 2.5 K and Hq = 0; 



moreover, the data are normalized to the maximum varia- 
tion, Z\i?™ aa: = Rn — R res . The continuous lines reported 
in the figures are the best-fit curves obtained by the model 
reported in Sec. |U 

In both samples, R s does not show any variation as 
long as the magnetic field reaches a certain value, depend- 
ing on T, that identifies the first-penetration field, H p . For 
Hq > Hp, vortices start to penetrate the sample and, con- 
sequently, R s increases. 

Figure [3] refers to the results obtained at T = 4.2 K. 
At this temperature, in both samples the R S (H ) curves 
exhibit a magnetic hysteresis, which disappears for Hq 
higher than a certain value, indicated in the figure as H'. 
The inset in panel (b) shows a minor hysteresis loop ob- 
tained by sweeping Hq from to 0.25 T and back. The 
field-induced variations of R s in sample P0 show some 
anomalies. Firstly, the application of a magnetic field of 
~ 1 T, which is about H C 2/15, causes a R s variation of 
~ 35% of the maximum variation. These field-induced 
variations of R s are much greater than those expected 
from the models reported in the literature 30,31,32 and 
detected in other superconductors 28,29]. A comparison 
with the results of panel (b) shows that in sample P6 a R s 
variation of the same order is obtained for the same value 
of H , even though, in this case, 1 T is about H c2 /2. Re- 
sults similar to those obtained in P0 have been observed in 
other MgB2 samples, produced by different methods and, 
therefore, seems to be a peculiarity of MgB 2 35,37,38, 
I40j . The finding that in sample P6 we have not observed 
this anomalous result strongly suggests that the enhanced 
R s variation is due to the two-gap superconductivity. 

A magnetic hysteresis of R s is expected in supercon- 
ducting samples in the critical state; it is ascribable to the 
different magnetic induction at increasing and decreasing 
DC fields |32]. Most likely, the different amplitude of the 
hysteresis loop obtained in the two samples is due to the 
different values of the critical current density; a smaller 
hysteresis is observed in sample P6 because of the smaller 
J c value. However, as it is visible in Fig. [3] also the shape 
of the hysteresis loop is different in the two samples. The 
decreasing- field branch of the R s (Hq) curve in sample P6 
has a negative concavity down to H p , as expected [42 . 
On the contrary, in sample P0 one can observe a plateau, 
in the field range 0.2 T, which cannot be justified in 
the framework of the critical-state models, considering the 
measured field dependence of J c (T3J . We would like to re- 
mark that this result has been obtained in all of the MgB2 
samples we have investigated [38|l40] . 

Figure |4] shows the field- induced variations of R s , for 
sample P0 (a), at T = 30 K, and for sample P6 (b), at T = 
7 K; for both samples, T/T c w 0.77. In the R S (H ) curve 
of sample P0 the hysteresis is still present, probably due 
to the high value of J c at this temperature, and still has an 
anomalous shape. We remark that in sample P0 we have 
observed magnetic hysteresis of R s up to T/T c w 0.95, 
while in sample P6 the hysteresis becomes undetectable 
at T/T c > 0.55. 

Figure [5] shows the field-induced variations of R s at 
temperatures near T c , where in both samples the R s (Hq) 
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»o H o ( T ) 

Fig. 3. Field-induced variations of R a for samples P0 (a) and 
P6 (b), at T = 4.2 K. AR S (H ) = R S (H ,T)-Rr ea , where R res 
is the residual mw surface resistance at T = 2.5 K and Hq = 0; 
AR™ ax = Rn — Rres- The line is the best-fit curve obtained, 
as explained in Sec. [5] with fioH C 2 = 1.71 T, uio/u = 0.67 and 
the field dependence of the critical current density reported in 
Ref. |13j . The inset shows a minor hysteresis loop obtained by 
sweeping Ho from to 0.25 T and back. 

curve is reversible. As will be shown in Sec. El at tem- 
peratures near T c the field dependence of the mw surface 
resistance can be accounted for by standard models also 
for sample P0, provided that the anisotropy of the upper 
critical field is taken into account. 

From isothermal R s (Hq) curves, obtained at different 
temperatures, we have deduced the temperature depen- 
dence of the characteristic fields, H p , H C 2 and H' . In Fig. [6] 
we report the values of H p , H C 2 and H' as a function of 
the reduced temperature, T/T° nset , for the two samples. 
The inset in panel (b) shows H C 2(T) of sample P6 in an 
enlarged scale. 

From Fig. [6^, one can see that H p of sample P0 ex- 
hibits a linear temperature dependence down to low tem- 
peratures, consistently with results reported by different 
authors in bulk [43|[44] and crystalline [45|[46] MgB 2 sam- 
ples. The extrapolated value at T = is about 55 mT; 
so, considering the demagnetization effect, the estimated 
value of the lower critical field is H c i(0) ss 70 mT. This 
value, although consistent with the lower critical field re- 
ported for MgB2 crystals by some authors [46], is slightly 



-j— i— i— i— |— i— t— t— r-i— i— i— r-T— i— i— |—t— t— t— r 

: Sample P0 




Vo ( T ) 

Fig. 4. Normalized field-induced variations of R 3 for samples 
P0 (a) and P6 (b), at T/T c ~ 0.77. The line in panel (b) is the 
best-fit curve of the data, obtained as described in the text by 
using the field dependence of the depinning frequency reported 
in Fig. H 

larger than that reported for bulk samples, which ranges 
from 15 to 45 mT. We suggest that this is ascribable to 
weak surface-barrier effects. 

In sample P6 we have obtained H p values smaller, but 
of the same order, than those of sample P0; this may 
be due to the irradiation effects. Indeed, the authors of 
Ref. [18] , from magnetization measurements in neutron- 
irradiated MgB2 crystals, have observed that the lower 
critical field reduces monotonically on increasing the flu- 
ence. 

The values of H c2 (T) indicated in Fig. [(J) as open tri- 
angles have been deduced measuring the magnetic field at 
which R s reaches the normal state value, R n . At the tem- 
peratures in which the upper critical field of sample P6 is 
higher than the maximum magnetic field achievable with 
our experimental apparatus (~ 1 T), the H C 2(T) values 
(full triangles in the figure) have been obtained as best- 
fit parameters, using the model reported in Sec. [4] On the 
contrary, since the results obtained in sample P0 cannot be 
accounted for by the model (except at temperature close 
to T c ), we report only the values we have directly deduced 
from the experimental results. For T > 5 K, the values we 
obtained for iJ c2 (T), in both samples, agree with those re- 
ported in Ref. p2] (the authors do not report H C 2 at lower 



6 



M. Bonura et al.: Microwave surface resistance in MgB2 



1.0 - 



0.9 - 



0.8 - 



51 



0.7 - 



t — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 



Sample P0 

7=38K n0$ 8j^$P O o 




p — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

Sample P6 

of 10 hr=8.5K 



0.9 - 



0.8 - 



0.7 - 
i i i 




(b)- 



50 100 150 200 



250 



Fig. 5. Normalized field-induced variations of R s for samples 
P0 (a) and P6 (b), at temperatures close to T c . Lines are best- 
fit curves, obtained by the model of Sec. [4] considering fluxons 
move in the flux-flow regime. The line of panel (a) has been 
obtained taking into account the anisotropy of the upper criti- 
cal field, as described in the text, with 7 = 3.3; for sample P0, 
we have set 7 = 1 consistently with the results of Ref. [14] , 



temperatures). Our results give complementary informa- 
tion about the temperature dependence of H C 2 of sample 
P6 at low temperatures. The continuous line in Fig. [SJd 
has been obtained by fitting the data of sample P6 with 
H c2 {T) = H c20 [l - (T/T c ) a }; we have obtained, as best-fit 
parameters, Mo#c20 = (2.2 ± 0.2) T, a = 1.9 ± 0.3 and 
T c = (8.9 ± 0.2) K. This temperature dependence of the 
upper critical field is consistent with that expected in con- 
ventional superconductors. On the contrary, for sample P0 
we observed an upward curvature of H c2 (T) , clearly visible 
in the inset, characteristic of two-gap MgB2 materials [TJ 

sua. 

H'(T) of Fig. [6J: corresponds to the value of the DC 
magnetic field at which the decreasing-field branch of the 
R s (Hq) curves deviates from the increasing- field branch. 
The zero values (without error bar) mean that the hystere- 
sis is not detectable at the corresponding temperatures. 
Consistently with the lower value of the critical current 
density, H'(T) is smaller in sample P6 than in P0. We 
would like to remark that the values of H'(T) could differ 
from the irreversibility field deduced from magnetization 
measurements. Indeed, it has been shown that, in sam- 
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Fig. 6. Temperature dependence of the characteristic fields, 
Hp, H C 2 and H', for the two samples. In panel (b): the inset 
shows the H C 2(T) values of sample P6 in an enlarged scale; the 
continuous line is the best-fit curve of the experimental data, 
obtained for sample P6, as described in the text. 



pies of finite dimensions, the application of an AC mag- 
netic field normal to the DC field can induce the fluxon 
lattice to relax toward an uniform distribution [47] . Fur- 
thermore, measurements we have performed in different 
superconducting samples have pointed out that, for sam- 
ples of millimetric size, the sensitivity of our experimental 
apparatus does not allow resolving magnetic hysteresis of 
R s when J c < 10 4 A/cm 2 . By considering the values of 
J c reported in Ref. [13], at T = 4.2 K we should obtain 
HqH' « 0.2 T for sample P6 and /iqH' w 3.5 T for sam- 
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pie PO. From Fig. [6b, one can see that this expectation is 
verified in sample P6; on the contrary, the value of H' in 
sample PO is about one order of magnitude smaller than 
the expected one. 



4 The model 

Microwave losses induced by static magnetic fields have 
been investigated by several authors 27,28,29,30,31,32. 
H2TH5] . At low temperatures and for applied magnetic fields 
lower enough than the upper critical field, the main con- 
tribution arises from the fluxon motion; however, it has 
been pointed out that a noticeable contribution can arise 
from the presence of normal fluid, especially at temper- 
atures near T c and for magnetic fields of the same order 
of H C 2(T). The majority of the models assume an uni- 
form distribution of fluxons inside the sample; so, they 
disregard the effects of the critical state. Very recently, 
we have investigated the field-induced variations of the 
mw surface resistance in superconductors in the critical 
state [12148] . and have accounted for the magnetic hys- 
teresis in the R S (H ) curves. 

In the London local limit, the surface resistance is pro- 
portional to the imaginary part of the complex penetra- 
tion depth, A, of the em field: 



R s = —flow Im[A(w, B,T)]. 



(1) 



The complex penetration depth has been calculated in 
different approximations [3"Ull31| . Coffey and Clem (CC) 
have elaborated a comprehensive theory for the electro- 
magnetic response of superconductors in the mixed state, 
in the framework of the two- fluid model of superconductiv- 
ity [30]. The theory has been developed under two basic 
assumptions: i) inter- vortex spacing much less than the 
field penetration depth; ii) uniform vortex distribution in 
the sample. With these assumptions vortices generate a 
magnetic induction field, B, uniform in the sample. This 
approximation is valid for Hq > 2H c i whenever the fluxon 
distribution can be considered uniform within the AC pen- 
etration depth. 

In the linear approximation, <C Hq, X(uj, B ,T) 
expected from the CC model is given by 



with 



\{B,T) 
8 nf {u,B,T) = 



An 



~ ^l-woni-BIB^T)]' 

5q(u) 

y/l-[l-w (T)]\l-B/B c2 (T)] 



(3) 
. (4) 



where Ao is the London penetration depth at T = 0, do 
is the normal-fluid skin depth at T = T c , Wq(T) is the 
fraction of normal electrons at Hq — 0; in the Gorter and 
Casimir two-fluid model w (T) = (T/T c ) 4 . 



S v is the effective complex skin depth arising from the vor- 
tex motion; it depends on the relative magnitude of the 
viscous and restoring-pinning forces, which identities the 
depinning frequency loq. 5 v can be written as 



<J2 Si V uj J 



where 



250c 



(5) 



(6) 



fiQLUT] 

with r\ the viscous-drag coefficient and 4>q the quantum of 
flux. 

When the frequency of the em wave, to, is much lower than 
luq, the fluxon motion is ruled by the restoring-pinning 
force. On the contrary, for w > woi the fluxon motion 
takes place around the minimum of the pinning-potential 
well and, consequently, the restoring-pinning force is nearly 
ineffective. So, the contribution of the viscous-drag force 
predominates and the induced em current makes fluxons 
move in the flux-flow regime. In this case, enhanced field- 
induced energy losses are expected. 

As it is clear from Eqs. (1-4), it is expected that the 
features of the R S (H ) curves strongly depend on the 
applied-field dependence of B. On the other hand, the CC 
theory is strictly valid when B is uniform inside the sam- 
ple; in particular for H ^> H&, the R S (H ) curves can be 
described setting B = ^oHq. When fluxons are in the crit- 
ical state, the assumption of uniform B is no longer valid 
and the CC theory does not correctly describe the field- 
induced variations of R s . As a consequence, the hysteresis 
in the R s (Hq) curve cannot be justified by Eqs. (1-4). In 
our field geometry (see Fig. [IJ,), the effects of the non- 
uniform B distribution on R s are particularly enhanced 
because in the two surfaces of the sample normal to the 
external magnetic field the mw current and fields pen- 
etrate along the fluxon axis and, consequently, the mw 
losses involve the whole vortex lattice. However, in this 
case, one can easily take into account the non-uniform B 
distribution by calculating a proper averaged value of R s 
over the whole sample as follows [421148] 



R. = Tg I HA Bir) )<IS. 



(7) 



where S is the sample surface, S is its area and r identi- 
fies the surface element. 

The pinning effects are particularly enhanced at temper- 
atures smaller enough than T c , where the dissipations are 
essentially due to vortex motion. So, the main contribu- 
tion to R s comes from the sample regions in which fluxons 
experience the Lorentz force due to the mw current, i.e. 
where Hq x J u ^ 0. Furthermore, in order to take into due 
account the critical-state effects by Eq. (J7J), it is essential 
to know the B profile inside the sample, determined by 
J C (B). 

Recently, using this method, we have investigated the 
effects of the critical state on the field-induced variation 
of R s , at increasing and decreasing fields [42[l48], We have 
shown that the parameter that mainly determines the pe- 
culiarities of the R s (Hq) curve is the full penetration field, 
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H*. Firstly, the width of the hysteresis is directly related 
to the value of H*; samples of small size and/or small J c 
are expected to exhibit weak hysteretic behavior. Further- 
more, H* determines the shape of the hysteresis loop as 
well. On increasing the external field from zero up to H* , 
more and more sample regions contribute to the mw losses; 
this gives rise to a positive curvature of the increasing- field 
branch of the R s (Hq-) curve. For Hq > H* , in the whole 
sample the local magnetic induction depends about lin- 
early on the external field and the increasing-field branch 
is expected to have a negative concavity. The shape of the 
decreasing-field branch is strictly related to the shape of 
the magnetization curve; it should exhibit a negative con- 
cavity, with a monotonic reduction of R s in the whole field 
range swept. 
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5 Discussion 

As we have shown in Sec. the R S (H ,T) curves exhibit 
different peculiarities in the unirradiated sample (P0) and 
the strongly irradiated sample (P6). The model described 
in Sec. [4] fully justifies the experimental results obtained 
in sample P6, which exhibits a single-gap superconductiv- 
ity [T5lfT6] . On the contrary, the results obtained in sample 
P0 cannot be justified in the framework of the same model, 
either for the magnetic-field dependence or for the tem- 
perature dependence of the surface resistance, even at zero 
DC field. Only the results obtained at temperatures close 
to T c can be justified, provided that the anisotropy of the 
upper critical field is taken into due account. In the fol- 
lowing, firstly we will discuss the temperature dependence 
of the mw surface resistance in the absence of DC mag- 
netic fields; successively, we will discuss the field-induced 
variations of R s . 

5.1 Temperature dependence of R s in zero magnetic 
field 

Figure [7] shows the normalized values of the surface resis- 
tance at Hq = as a function of the reduced temperature 
for both samples. The R S (T) curve of sample P6 shows 
a wide transition, broadened in a roughly symmetric way 
with respect to the middle point at R s /R n = 0.5. This 
behavior can be ascribed to the T c distribution over the 
sample. On the contrary, in sample P0 one can notice a 
sharp variation of R S (T), at temperatures near T c , and a 
wide tail, extending from T/T c w 0.9 down to T/T c « 0.7, 
which cannot be ascribed to the T c distribution. The lines 
in the figure are best-fit curves; they have been obtained 
with different procedures for the two samples. 

The results obtained from the model discussed in Sec. 2] 
setting B — in Eqs. (HHi]) converge to those of the em 
response of superconductors in the Meissner state, in the 
framework of the two-fluid model. In this case, the tem- 
perature dependence of R s /R n is determined, apart from 
the T c distribution over the sample, by the temperature 
dependence of the normal- fluid density, wq(T), and the ra- 
tio Xq/Sq- In order to fit the experimental data obtained 



Fig. 7. Normalized values of the mw surface resistance as a 
function of the reduced temperature, obtained in the two sam- 
ples at Ho = 0. Symbols are experimental data; lines are the 
best-fit curves obtained as described in the text. 

in sample P6, we have assumed wq(T) = (T/Tc) 41 , consis- 
tently with the Gorter and Casimir two-fluid model, and 
have used Eqs. fTHI|) with B — 0. We have averaged the 
expected curve over a gaussian distribution function of T c 
with T c0 = 8.5 K and a = 0.2 K (see Sec. [2]) and have used 
Ao/<5o as fitting parameter. The best-fit curve, dashed line 
in Fig. [3 has been obtained with Ao/^o = 0.14; however, 
we have found that the expected curve is little sensitive 
to variations of Ao/^Oj except at low temperatures, where 
the measured R 8 is limited by the sensitivity of our ex- 
perimental apparatus. In particular, by varying T c q and 
a within the experimental uncertainty, good agreement is 
obtained with Xo/5q values ranging from 0.04 to 0.15. This 
occurs because the T c distribution broadens the R S {T) 
curve, hiding the Xq/Sq effects. 

Unlike for sample P6, the results of Fig. [7] obtained 
in P0 cannot be justified in the framework of the Gorter 
and Casimir two-fluid model, using reasonable values of 
Xo/5o- On the other hand, different authors [2115114^] have 
shown that the temperature dependence of the field pen- 
etration depth in MgB2 cannot be accounted for by ei- 
ther the Gorter and Casimir two-fluid model or the stan- 
dard BCS theory. A linear temperature dependence of the 
condensed fluid density, in a wide range of temperatures 
below T c , has been reported, which has been justified in 
the framework of two-gap models for the MgB2 supercon- 
ductor [51l4"9"]. Prompted by these considerations, we have 
hypothesized a linear temperature dependence of wq. The 
continuous line in Fig. [7] is the best-fit curve; it has been 
obtained by Eqs. (fTHU) with B = 0, w (T) = T/T c and 
Xo/So — 0.15. The wide low-T tail is essentially deter- 
mined by the linear temperature dependence of wq. The 
sensitivity achievable by our experimental apparatus does 
not allow determining the small variations of R S {T) for 
T/T c < 0.5; so, from R S (T) measurements no indication 
about the temperature dependence of the densities of the 
normal and condensed fluids at low temperatures can be 
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obtained. However, the linear temperature dependence of 
the lower critical field we obtained (see Fig. [6^) strongly 
suggests that wq linearly depends on T down to low tem- 
peratures. 

5.2 Field dependence of R s in sample P6 

In conventional (single-gap) superconductors, it is expected 
that the field dependence of the mw surface impedance is 
described by the model reported in Sec. HI In this frame- 
work, in order to calculate the expected field-induced vari- 
ations, by Eqs. JIHI]), the essential parameters are the 
value of Xq/Sq, H C 2(T), the depinning frequency, luq, and 
its field dependence. It is not necessary to consider the 
upper-critical-field anisotropy, 7, because it has been shown 
that in the P6 sample 7 = 1 [33] . When the critical-state 
effects cannot be neglected, in order to use Eq. (J7]), it is 
also essential to know the profile of the induction field 
determined by the field dependence of the critical current 
density. The value of Ao /So has been determined by fitting 
the R S (T) curve at H = 0; the critical current density and 
its field dependence are reported in Ref. [T3j. The values 
of H c2 (T) at T > 5 K are reported in Ref. [13], and/or 
deduced from our experimental data; at T < 5 K, H c2 has 
to be considered as parameter. It is worth noting that the 
large uncertainty of Ao/<5o, we obtained for this sample, 
does not affect the best-fit curves because this parameter 
essentially determines the normalized AR s (Hq) value at 
B = 0. 

As one can see from Fig. [31 at T = 4.2 K the R S {H ) 
curve exhibits a magnetic hysteresis, indicating that the 
effects of the critical state are not negligible. In order to 
use Eq. ([7]) , we have calculated the B profile in the sample 
using the field dependence of the critical current, J C (B), 
reported in Ref. [13j and we have set the induction field at 
the edges of the sample as B = /j,q(Ho—H p ); we have taken 
H C 2 and uio as fitting parameters. The line of Fig. [3b is the 
best-fit curve; it has been obtained with /j,qH C 2 — 1.71 T 
and ujq/ui = 0.67 independent of Ho- For the sake of clar- 
ity, in Fig. [8]we report the results obtained by sweeping the 
DC magnetic field from to 0.25 T and back, along with 
the expected curve. The inset shows the B profile along the 
width of the sample at half height, determined by J C (B); 
the continuous lines are the increasing-field profiles, the 
dashed ones are the decreasing-field profiles at the same 
external-field values. As one can see, taking into account 
the field distribution inside the sample, the experimen- 
tal results are quite well justified in the framework of the 
model discussed in Sec.[U In the increasing-field branch, a 
change of concavity is well visible at fj,o(Ho—H p ) ss 0.04 T, 
consistently with the expected value of the full penetra- 
tion field; the decreasing-field branch exhibits a negative 
concavity in the whole range of fields. 

Following the procedure above described, from the best 
fit of the experimental data of the isothermal R s (Ho) 
curves at T < 6 K we have obtained the H C 2(T) values 
indicated as full triangles in Fig. [6Jd. 

When the R S (H ) curves do not show hysteresis, the 
effects of the critical state are negligible and the induction 




Fig. 8. Field-induced variations of the mw resistance, obtained 
in sample P6 by sweeping the magnetic field from to 0.25 T 
and back. The line is the best-fit curve obtained, as explained 
in the text, with /io-ff C 2 = 1.71 T, u>q/u = 0.67 and the field 
dependence of the critical current density reported in Ref. [13j . 
The inset shows the B profile at increasing ( — ) and decreasing 
(- -) fields; w is the width of the sample. 

field, B, can be considered uniform. In this case, we have 
used the following approximate expression for the magne- 
tization: 

M = -H p + Hp (Ho - H p ) ; 

tl c 2 — tip 

and, consequently 

B = Mo (l + H f* Hp ) ( H o - h p)- 

Several calculations have shown that, in order to fit the 
experimental data at T > 7 K, it is essential to consider 
the T c distribution over the sample. So, we have aver- 
aged the expected R s (Hq) curves [calculated by Eqs. (JTJ- 
HJ] over the T c distribution (see Seed]). We have used for 
H P (T) and H C 2(T) the values deduced from the experi- 
mental data, letting them vary within the experimental 
uncertainty, and have considered the depinning frequency 
as parameter. The lines of Figs. H] and [5] have been ob- 
tained by this procedure. 

By fitting the results at T = 7 K, we have obtained the 
field dependence of wo/w reported in Fig. (5) The roughly 
constant value of ujo/lo we obtained up to [IqHo ~ 0.3 T 
indicates that in this field range individual vortex pinning 
occurs; on further increasing the magnetic field, the in- 
teraction between fluxons becomes important, collective 
vortex pinning sets in and, consequently, the depinning 
frequency decreases. The data obtained at hoHq > 0.65 T 
are well fitted setting iool^ — in Eq. ([5]); this means 
that, at T = 7 K and H0H0 > 0.65 T, the induced mw 
current makes fluxons move in the flux-flow regime. 

The best-fit curve of Fig. [5] has been obtained with 
wq/w = 0, as expected. Indeed, at temperature very close 
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Fig. 9. Magnetic field dependence of the depinning frequency, 
obtained for sample P6 by fitting the experimental results re- 
ported in Fig. [4b. 



to T c , the pinning effects are weak and the induced mw 
current makes fiuxons move in the flux-flow regime. 



5.3 Field dependence of R s in sample P0 

It has been shown by several authors that the proper- 
ties of the two-gap MgB2 superconductor in the mixed 
state cannot be accounted for by standard theories [T][31 
18,19,20,35,36,39]. It is by now accepted that this is re- 
lated to the double-gap nature of MgB 2 that is responsi- 
ble for an unusual vortex structure. Indeed, it has been 
highlighted, both experimentally and theoretically, that 
the vortex cores are characterized by two different spatial 
and magnetic-field scales [2T1I22] . Because of the different 
magnetic-field dependence of the two gaps, on varying the 
field, the structure of the vortex lattice changes in an un- 
usual way. At low magnetic fields, quasiparticles by 7r and 
a bands are trapped within the vortex core, even if on 
different spatial scales because of the different coherence 
lengths £„• and in the field range 0.5 1 T, though a- 
band quasiparticles remain localized, the 7r-band quasipar- 
ticles spread over the sample |22) ; on further increasing the 
field, A n is strongly reduced, the 7r-quasiparticle contribu- 
tion remains almost unchanged while the cr-quasiparticlc 
contribution continues to increase with about the same 
rate up to the macroscopic H c i- So, a further charac- 
teristic field is needed for determining the fluxon-lattice 
properties of MgB 2 ; the existence of this crossover field, 
often indicated as H^ 2 , has been highlighted in several 
experiments 1,3, 22 , 23 , 24, 25 . The field-induced evolu- 
tion of the vortex lattice is expected to affect the vortex- 
vortex and vortex-pinning interactions, making the stan- 
dard models most likely inadequate to describe the fluxon 
dynamics. 

Results on the field-induced variations of R s in MgB2 
have been reported by some authors [5M5fjll3Tl l38 , 39 . 40 . 
[50] : most of them have highlighted several anomalies, which 



cannot be explained in the framework of standard mod- 
els for fluxon dynamics. In particular, it has been high- 
lighted unusually enhanced field-induced mw losses at ap- 
plied magnetic field much lower than H C 2- Only Zaitsev et 
al. [50] have explained the frequency and field dependence 
of the mw surface resistance of MgB2 films in the frame- 
work of standard models. The results we have obtained 
in sample P0 are similar to those reported by Shibata et 
al. [35] , who investigated the field dependence of the sur- 
face impedance in MgB2 single crystal in a wide range of 
DC magnetic fields (up to 14 T). At low temperatures, 
the authors have observed an initial fast variation of the 
field-induced mw dissipation up to fields of the order of 
1 T, followed by a slower one at higher fields. Consis- 
tently with the sharp field-induced variation of the heat 
capacity and thermal conductivity, the enhanced low-field 
variation of the mw losses has been ascribed to the high 
increase of ir quasiparticles in the vortex cores. At higher 
fields, the variation is slower because of the saturation of 
the 7r-quasiparticle contribution. 

According to Shibata et al., the enhanced field-induced 
variation we observed in sample P0 can be qualitatively 
ascribed to the strong reduction of A v in the field range we 
have investigated. However, the observed R s (Hq) curves 
differ from the expected ones in both the intensity and the 
shape. Here, we discuss the shape of the R S (H ) curves. 

In a wide range of temperatures below T c , we have ob- 
served a magnetic hysteresis, which should be related to 
the different magnetic induction at increasing and decreas- 
ing fields, due to the critical state. As discussed in Sec. [4] 
the increasing- field branch of the R s (Hq) curve should 
exhibit a change of concavity, from positive to negative, 
when the external magnetic field reaches the full pene- 
tration field, H*. By considering the sample width and 
the value of J c at T = 5 K reported for sample P0 [TB] . 
the expected value of H* is « 2.6 T. Nevertheless, we 
observe a negative concavity of the R s (Hq) increasing- 
field branch in the whole range of fields investigated (see 
Fig. [3^) , even if the maximum value of the applied field is 
well below H*. The decreasing- field branch should show 
a monotonic reduction of R s down to low fields. In con- 
trast, for Hq < H' , we observe an initial weak reduction 
followed by a plateau, from /io-ffo ~ 0-2 T down to zero. 
The presence of this plateau is puzzling because it would 
suggest that the trapped flux does not change anymore on 
decreasing the field below ~ 0.2 T, although this value is 
four times larger that H p . 

Another anomalous result concerns the range of mag- 
netic fields in which we observe the hysteretic behavior. 
As we have already mentioned, we have experienced that 
for samples of millimetric size the sensitivity of our exper- 
imental apparatus allows detecting hysteresis in R S (H ) 
f° r Jc ^ 10 4 A/cm 2 . From Fig. 9 of Ref. [13], one can de- 
duce that, in sample P0, such condition occurs at ^o^o ~ 
4 T; so, we should detect hysteresis in the whole range of 
fields we have investigated. On the contrary, we obtained 
if' (4.2 K) ~ 0.5 T, one order of magnitude lower than the 
expected value. 
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We would like to remark that these anomalies have 
been observed in all the bulk MgB2 samples (unirradiated) 
we have investigated, no matter the preparation method 
and the components ( 11 B or 10 B) used in the synthesis 
process |381l40j . The finding that in sample P6 the ex- 
perimental results are fully justified by the used model, 
strongly suggests that these anomalies are strictly related 
to the presence of the two superconducting gaps. 

At temperatures close to T c and for H > 0.5H C 2(T), 
the experimental results can be accounted for by the model 
discussed in Sec. HI provided that the anisotropy of the 
upper critical field is taken into due account. Following 
Ref. [14] . to take into account the anisotropy, we have as- 
sumed that the polycrystalline sample is constituted by 
grains with the c-axis randomly oriented with respect to 
the DC-magnetic-field direction; so, the distribution of 
their orientations follows a sin(0) law, being 9 the angle 
between Hq and c. Furthermore, we have used for the an- 
gular dependence of the upper critical field the anisotropic 
Ginzburg-Landau relation 



H c2 (6) 



sff cos 2 (0) + sin 2 (0) 



where 7 = H^f/H^ is the anisotropy factor. 

The field-induced variations of R s observed at T = 
38 K (reported in Fig.[5k) do not exhibit hysteresis; so, in 
this case, B can be considered uniform. Furthermore, at 
temperatures near T c , one can reasonably suppose fluxons 
move in the flux-flow regime. In this condition, the ex- 
pected R S (H , H c2 (9)) curve depends on Ao/<5o (obtained 
by fitting the R S (T) curve at H = 0), H^ c and 7. On 
the other hand, the H c2 values deduced from the isother- 
mal R S (H ) curves (see Fig. EJd) coincide with the mag- 
netic field at which the whole sample goes to the normal 
state, i.e. H^f. In order to fit the results at T — 38 K, we 
have averaged the expected curve [calculated by Eqs. ((Tr- 
ill)] over a sin(6 l ) distribution, have used for H^f the value 
of the magnetic field at which R s /R n — 1, letting it vary 
within the experimental uncertainty, have taken 7 as free 
parameter. At this value of temperature, the experimen- 
tal results can be accounted for using 7 = 3.3±0.5. In 
particular, the best-fit curve reported in Fig. [5^ has been 
obtained with 7 = 3.3 and n^H^f = 145 mT. As one can 
see, the field-induced variation of R s is well described by 
the model of Sec. |H We think that this occurs because 
at this temperature the superfluid fraction of the tt band 
is strongly suppressed at low magnetic fields, the flux line 
gets a conventional structure and the fluxon dynamics can 
be described by standard models. 

Prompted by the results obtained at T — 38 K, we 
have tried to fit the experimental data obtained in the 
temperature range 34 37 K by the same method (for 
these temperatures, the upper critical field has been di- 
rectly deduced from the R S (H ) curves). Since in this 
temperature range we have detected magnetic hystere- 
sis, we have considered only the reversible part of the 
R s (Hq) curve. In order to fit the data, we have hypothe- 
sized fluxons move in the flux-flow regime, have considered 
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Fig. 10. Normalized field-induced variations of R s for sam- 
ple P0, at different temperatures near T c . Open symbols are 
the results obtained at increasing Ho, full symbols those at 
decreasing Ho- The R s (Ho) curve at T — 38 K is reversible. 
The lines are the expected curves, obtained with 7 = 2.6 as 
described in the text. The inset shows the temperature depen- 
dence of the magnetic field H cr , above which the experimental 
data can be justified in the framework of the model of Sec. [4] 



for H^f(T) the values reported in the inset of Fig.[6jD, let- 
ting them vary within the experimental uncertainty, and 
have taken 7 as fitting parameter. We have found that 
at high fields the experimental results can be fitted us- 
ing 7 = 2.6 ± 0.2. Fig. [10] shows a comparison between 
the expected curves, obtained with 7 = 2.6, and the ex- 
perimental data for T = 34 38 K; open symbols are the 
results obtained at increasing Hq, full symbols those at de- 
creasing Hq. As one can see, the expected R S (H ) curve at 
T = 38 K so obtained poorly agrees with the experimen- 
tal data at low fields; on the contrary, the line of Fig. [3^, 
which has been obtained with 7 = 3.3, fits the data in the 
whole range of magnetic fields. However, the value 7 = 2.6 
is closer to the upper-critical-field anisotropy reported in 
the literature for MgB2 at temperatures near T c 



The results of Fig. [TTJ show that for Hq greater than 
a certain threshold value, depending on T , the data can 
be justified in the framework of the model describing the 
fluxon dynamics of conventional vortex lattice. The tern- 
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perature dependence of the threshold field, H cr , is re- 
ported in the inset. For Hq < H cr the field-induced mw 
losses are larger than those expected for single-gap su- 
perconductors in the mixed state. We suggest that this 
surplus of mw losses is due to the additional contribu- 
tion of 7r-band quasiparticles within the vortex cores with 
respect to that of one-gap superconductors; the finding 
that H cr decreases on increasing T seems to support this 
hypothesis. It is easy to see that the H cr (T) values coin- 
cide, within the experimental uncertainty, with 0.5H C 2(T) 
(see the inset of Fig. [5]). Presently, it is not clear why 
just above 0.5iJ C 2 the results can be justified by a stan- 
dard model, which does not consider the two-gap nature of 
MgPj2. Furthermore, we remark that H cr cannot be iden- 
tified with the magnetic field at which the 7r-band super- 
fluid is suppressed; indeed several authors have reported 

6 Conclusions 

We have investigated the microwave surface resistance at 
9.6 GHz of two polycrystalline Mg 11 B2 samples prepared 
by direct synthesis from Mg (99.999% purity) and crys- 
talline isotopically enriched n B (99.95% purity). That la- 
belled as PO consists of pristine material; the other, la- 
belled as P6, has been exposed to neutron irradiation 
at very high fluence. Several superconducting properties 
of these samples have been reported in Refs. [T3IIT4. 15. 
[T6lll7j. Point-contact spectroscopy and specific- heat mea- 
surements, have shown that sample PO exhibits a clear 
two-gap-supercon-ductivity behavior; in sample P6 the ir- 
radiation process determined a merging of the two gaps 
into a single value. To our knowledge, the mw response of 
neutron irradiated MgPj2 samples has not yet been inves- 
tigated. 

The mw surface resistance has been measured as a 
function of the temperature and the DC magnetic field. By 
measuring the field-induced variations of R s at increasing 
and decreasing fields we have detected a magnetic hystere- 
sis ascribable to the critical state of the fluxons lattice. 
The range of temperatures in which the hysteretic behav- 
ior has been observed is different for the two samples; 
in the irradiated sample the hysteresis is undetectable at 
T/T c > 0.55 while in the unirradiated sample it is de- 
tectable up to T/T c w 0.95. 

The results obtained in the irradiated sample have 
been quite well justified in the framework of the Coffey 
and Clem model with the normal fluid density following 
the Gorter and Casimir two-fluid model. In order to ac- 
count for the hysteretic behavior, we have used a gen- 
eralized Coffey and Clem model in which we take into 
account the non-uniform fluxon distribution due to the 
critical state. 

The peculiarities of the mw surface resistance of sam- 
ple PO differ from those observed in sample P6, in both the 
temperature and the field dependencies. The R S {T) curve 
obtained at zero field shows a wide tail, from T/T c 0.9 
down to T/T c w 0.7, which cannot be justified in the 
framework of the Gorter and Casimir two-fluid model. We 



have shown that, in order to account for this behavior, it is 
essential to hypothesize a linear temperature dependence 
of the normal and condensed fluid densities. Such finding 
agrees with the experimental temperature dependence of 
the penetration depth reported in the literature, which 
have been justified in the framework of two-gap models 
for the MgB2 superconductor. 

The R s (Hq) curves in sample PO have shown several 
anomalies, especially at low temperatures, among which 
an enhanced field-induced variation and a magnetic hys- 
teresis of unconventional shape. At low temperatures, a 
magnetic field Hq H C 2/15 causes a R s variation of 
« 35% of the normal-state value. We remark that in sam- 
ple P6 a variation of the same order of magnitude is ob- 
tained for Hq k H C 2/2. The shape of the magnetic hys- 
teresis, which has been observed in a wide range of tem- 
peratures below T c , cannot be justified in the framework 
of the critical-state models; the most unexpected behavior 
concerns the decreasing-field branch, in which we observed 
a plateau extending from /io-ffo ~ 0.2 T down to zero. The 
presence of this plateau is puzzling because it would sug- 
gest that the trapped flux does not change anymore on 
decreasing the field below 0.2 T, although this value is 
four times larger than the first penetration field. 

The investigation at temperatures near T c has high- 
lighted that, in the range T = 34 -f- 38 K, the results 
obtained in sample P0 for Hq > 0.5H C 2 can be justified 
in the framework of the Coffey and Clem model taking 
into account the anisotropy of the upper critical field. We 
suggest that this occurs because at these field values the 
superfluid fraction of the 7r band is strongly suppressed, 
the flux line gets a conventional structure and the fluxon 
dynamics can be described by standard models. 

The enhanced field- induced variation of R s , observed 
at low T in the whole range of fields investigated as well as 
at T ~ T c for Hq < 0.5H C 2, may be qualitatively ascribed 
to the presence and motion of the giant cores due to the 
7r-band quasiparticles. On the contrary, the origin of the 
anomalous shape of the R s (Hq) curve is so far not un- 
derstood. We would like to remark that the results we ob- 
tained in sample P0 are very similar to those, not reported 
here, we have obtained in several MgB2 samples (unirra- 
diated), no matter the preparation method and the com- 
ponents ( U B or W B) used in the synthesis process. The 
comparison between the results obtained in the two sam- 
ples here investigated strongly suggest that the anomalies 
in the R s {Hq) curves are related to the unusual structure 
of fluxons due to the two superconducting gaps. According 
to what suggested by different authors, our results confirm 
that the standard models are inadequate to describe the 
fluxon dynamics in two-gap MgB2. Further investigation 
is necessary for understanding how to take into account 
the complex vortex structure in describing the fluxon dy- 
namics in MgB2. 
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